Bubble dynamics: (nucleating) radiation inside dust 



O 

(N 

G 



R. Casadio 1 * 2 -* and A. Orlandi 1 ' 2 ' t 

1 Dipartimento di Fisica, Universita di Bologna, via Irnerio 46, 4-0126 Bologna, Italy 
2 INFN, Sezione di Bologna, Via Irnerio 46, 1-40126 Bologna, Italy 

We consider two spatially flat Friedmann-Robertson- Walker spacetimes divided by a time-like 
thin shell in the nontrivial case in which the inner region of finite extension contains radiation and 
the outer region is filled with dust. We will then show that, while the evolution is determined 
by a large set of constraints, an analytical description for the evolution of the bubble radius can 
be obtained by formally expanding for short times after the shell attains its minimum size. In 
particular, we will find that a bubble of radiation, starting out with vanishing expansion speed, can 
be matched with an expanding dust exterior, but not with a collapsing dust exterior, regardless of 
the dust energy density. The former case can then be used to describe the nucleation of a bubble 
of radiation inside an expanding dust cloud, although the final configuration contains more energy 
than the initial dust, and the reverse process, with collapsing radiation transforming into collapsing 
dust, is therefore energetically favored. We however speculate a (small) decaying vacuum energy 
or cosmological constant inside dust could still trigger nucleation. Finally, our perturbative (yet 
analytical) approach can be easily adapted to different combinations of matter inside and outside 
the shell, as well as to more general surface density, of relevance for cosmology and studies of defect 
formation during phase transitions. 

PACS numbers: 04.40.-b,04.20.-q,04.20.Cv,98.80. Jk 
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I. INTRODUCTION 

The dynamics of an infinitely thin, spherically sym- 
metric shell S separating two spacetime regions Q± with 
given metrics is a well known problem of General Rela- 
tivity. The general theory dates back to 1965 [1] and is 
completely understood. Given the symmetry of the sys- 
tem, we can use spherical coordinates x± — {t±, r±, 9, <fi} 
in ft±, respectively, where the angular coordinates are 



the same in both patches, and < 



< r_ 



< r+, 



with r± — r±(t±) the (in general time-dependent) radial 
coordinates of the shell in Q,±. One then takes specific 
solutions g^ u of the Einstein equations inside f2± and im- 
poses suitable junction conditions across S, namely the 
metric is required to be continuous across the shell, 
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(1) 



whereas the extrinsic curvature K^j of £ is allowed to 
have a jump proportional to the surface stress-energy 
tensor of the time-like shell <7y (italic indices run on the 
shell's three-dimensional world-sheet), 

[K i j }-5 i j [K l l }=Ka i j , (2) 

in which [K^] = K i 3 \ + — K^]— denotes the difference 
between extrinsic curvatures on the two sides of the shell. 
Note that we set c = 1 and k = 8 7tGn/3 = £p/Mp (= 1 
when convenient), where Gn is Netwon's constant and £p 
(Mp) the Planck length (mass). Although the classical 
evolution equation (2) may appear simple, it has been 
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solved only in a few special cases, most notably for the 
vacuum or with cosmological constants in f2± [2] (for an 
extensive bibliography see Ref. [3]). 

A most intriguing result emerges in the semiclassical 
picture [4, 5], where one finds that "bubbles" can be 
quantum mechanically created from nothing (in a sense, 
at the expense of gravitational energy). This may oc- 
cur when one has a classical solution for an expanding 
shell areal radius with a (finite) minimum value (turn- 
ing point of the classical trajectory, larger than £p), and 
a non-vanishing quantum mechanical amplitude for the 
"tunneling" into such a system from one without the 
shell (that is, a shell of zero area). It has been conjec- 
tured that these bubbles could represent child universes 
generated inside a parent (or "landscape") universe [6— 
8] , if they expand indefinitely (or at least long enough) . 
Bubble dynamics might also be used to model regions 
of space within which a matter phase transition occurs 
(from false to true vacuum, as well as between different 
form of matter [9]). One can, for instance, use such a 
model to approximate the formation of radiation from a 
decaying scalar field during reheating after inflation. It 
is known that for an inflaton with a quadratic potential, 
the time averaged dynamics of the final oscillation phase 
mimics that of matter. The approach developed in this 
paper could turn out to be suitable to describe the decay 
into radiation. In the end, knowing the correct evolution 
of such a bubble would be of great help in understanding 
how defects formed during phase transitions are "ironed 
out" by the expansion of the new phase. 

In this paper, we are mainly interested in presenting 
an analytical (perturbative in time) approach to study a 
time-like shell's dynamics and to obtain analytical condi- 
tions for the existence of expanding bubbles in terms of 
the energy densities inside and outside the shell, when 



2 



such regions contain homogeneous dust or radiation. 
This problem is made technically cumbersome because 
of the occurrence of algebraic constraints (to ensure the 
arguments of proliferating square roots are positive) . We 
shall therefore consider in details only the specific case 
of nucleation of a spatially fiat radiation bubble inside 
(spatially flat) collapsing or expanding dust, in order to 
keep the presentation of our method more streamlined. 
Nonetheless, these cases are also of particular physical in- 
terest. For example, one can conceive the density inside a 
collapsing astrophysical object might be large enough to 
allow for the creation of supersymmetric matter which, in 
turn, would then annihilate regular matter and produce 
a ball of radiation [10]. Likewise, in a dense matter- 
dominated expanding universe, one might consider the 
possibility of spontaneous nucleation of radiation bub- 
bles. We shall find the dust energy density just sets 
the overall scale of the problem. Assuming the bubble 
surface density is (initially) constant, expanding radia- 
tion bubbles may then be matched with an expanding 
dust exterior, the time-like shell surface density being 
uniquely related to the inner radiation density. In order 
to view the bubble creation as a phase transition from 
dust to radiation (plus the surface density of the shell), 
one however needs an external source of energy, since the 
total energy of the bubble is larger than the initial en- 
ergy of the dust. This implies that the reverse process 
of collapsing radiation turning into collapsing dust is ac- 
tually favored energetically. Moreover, no configuration 
with expanding bubbles is allowed inside collapsing dust, 
regardless of how large is the dust energy density, and 
the conversion between the two types of matter therefore 
appears highly disfavored in this case. 

In Section II, we briefly review the fundamental equa- 
tions and constraints that describe general time-like bub- 
ble dynamics, following Ref. [9], and then specify all 
expressions for ± given by spatially flat Friedmann- 
Robertson- Walker (FRW) regions filled with dust or ra- 
diation. In Section III, we work out the explicit case 
of a radiation bubble of constant surface density nucle- 
ated inside collapsing or expanding dust, for which we 
obtain the initial minimum radius in terms of the inner 
and outer energy densities. Finally, in Section IV, we 
make some considerations about our findings and possi- 
ble future generalizations. 



II. BUBBLE DYNAMICS 

For our analysis, we will mostly follow the notation of 
Ref. [9], where the metric in each portion fi± of spacetime 
is given by 

d.s 2 = e"™ At 2 - e x ^ dr 2 - R 2 (r, t) dtt 2 , (3) 

where t = t± are the time coordinates inside the corre- 
sponding patches, and likewise for the three spatial co- 
ordinates. On the shell time-like surface E, one has the 



line element 

d S 2 | s = dT 2 -p 2 (r)df! 2 , (4) 

in which r is the proper time as measured by an observer 
at rest with the shell of areal radius p = R±(r±(t±),t±). 
The relation between r and t± is obtained from the equa- 
tion of continuity of the metric, Eq. (1), and is displayed 
below in Eq. (25) for the cases of interest. On solving 
Eq. (2) in terms of p and p = dp/ 'dr, one gets the dy- 
namical equation 

P 2 (t)=B 2 (t) P 2 (t)-1, (5) 

where 

(e + + 6_+9«q 2 /4) 2 -4e_ £+ 

= 9^ ' (6) 

with a — ct° (t) the shell's surface density and e± = 
e ±(t±) the time-dependent energy densities in tt± respec- 
tively. It is important to recall that metric junctions can 
involve different topologies for f2±, but we are here con- 
sidering only the so-called "black hole" type, in which 
both portions of spacetime have increasing area radii R± 
in the outward direction (of increasing r±). Assuming 
the surface density of the shell is positive, one must then 
have 1 

e + (r)-e_(r)> ^(r) , (7) 

at all times, in order to preserve the chosen spacetime 
topology [9, 11]. 

In the pure vacuum case, e± are constant and for con- 
stant a the solution is straightforwardly given by 

p(r) = B- 1 cosh(Sr) , (8) 

where B — B(e±,a) from Eq. (6) is also constant [2]. 

In the non-vacuum cases, finding a solution is how- 
ever significantly more involved. Regardless of the mat- 
ter content of Q±, it is nonetheless possible to derive a 
few general results for a bubble which nucleates at a time 
t = r , that is a shell that expands from rest, 

Po = , (9) 

with initial finite (turning) radius (po > 0), where the 
subscript will always indicate quantities evaluated at 
the time r = r . First of all, from Eq. (5), the initial 
radius must be given by 

Po = |S _1 | - (10) 
which requires B real, or 

(e + + £o- +9 kctq/ 4 ) 2 > 4e -eo+ • (H) 



1 This also implies that eo + > erj— and, from the Friedmann equa- 
tion (16) given below, > H'^. 
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This condition is always satisfied if e dust and eg ad are 
both positive and will therefore be of no relevance in this 
paper, but must be carefully considered when allowing 
for negative energy densities (and non-vanishing spatial 
curvature). Further, upon deriving Eq. (5) with respect 
to r (always denoted by a dot) 

2pp = 2(BBp 2 + B 2 pp^j , (12) 
and using Eq. (9), one also obtains 

Bo = , (13) 
assuming p is not singular. The constraint in Eq. (7) at 

T = TO, 



(14) 



and the conditions in Eqs. (9) and (13) will play a crucial 
role in the following. 



A. Flat FRW regions 

Since we wish to study the particular case of a shell E 
separating two regions Q± filled with homogeneous fluids, 
the metrics in Vt± will be taken to be the usual FRW 
expressions. Moreover, we already assumed f2_ has finite 
initial extension and, by definition, represents the interior 
of the shell. As a further simplification, we shall only 
consider flat spatial curvature and set the cosmological 
constant A = everywhere. 

The metrics (3) in the inner and outer regions are 
therefore given by 

d S 2 =dt 2 -a 2 (t)[dr 2 +r 2 (d9 2 +sm 2 ed ( f> 2 )] , (15) 

where a(t) is the scale factor which evolves according to 
the Friedmann equations 





/l da 


H 2 = 






\a dt 


n 1 d 2 a 


(I da 


2 ad^ +i 


{ a dt 



: K e 



= —3np 



(16) 
(17) 



We assume the energy density e and pressure p of the 
fluids obey barotropic equations of state, 

p = we(t), (18) 

and recover the well-known behaviors 



e(t) 



a(t) 
a 



3(w+l) 



= e 



(19) 



in which eo is the density evaluated at a reference instant 
of time t = to and a = a(t ). For dust, w = (p = 0), 



whereas for radiation w = 1/3, so that 

3 

,dustfj.\ _ £ fl 
( ) " 3(t) 

a*(t) ' 



(20) 



e rad (i) 



The evolution of scale factors in cosmic time for expand- 
ing (f) and contracting (|) solutions are finally given by 



a\f(t) = (7 ± 2 VM rad tj 



1/2 



da^ d y/M^d 



(21) 



dt 



a rad (t) 



and 



^dust 

2 n 



(t) = (s±% VW^ tj 



2/3 



da n st 

dt 



± 



I M dust 
a dust (t) ' 



(22) 



where, in the above r.h.s., the + signs are for expansion 
and — signs for contraction, 



M rad = itaj e£ ad 



M d 



ust 



-.3 ^dust 



(23) 



and 7 and 5 integration constants that determine the 
size of the scale factors at t = 0. Later, for convenience, 
we will set 7 = S = 1 at t = 0, so that e(0) = e and 
a(0) = ao = 1. 

Let us now consider the time-like shell S at r± = 
r ±(^±) separating the two regions Q±. Clearly, metric 
continuity implies 



p = a±(t±) r s ± (t ± ) . 



(24) 



The inner and outer spaces are characterized by different 
physical parameters. In particular, as one can see from 
Eq. (6), the shell's dynamics are determined by: 

1. The type of fluid inside the shell (its equation of 
state W-); 

2. The initial values ao_ and e _; 

3. The type of fluid outside the shell (its equation of 
state w + ): 

4. The initial values ao+ and eo+; 

5. The shell surface density a (as a function of the 
radius p). 

A given configuration of dust, radiation and surface den- 
sity is admissible only if the corresponding initial condi- 
tions are such that Eqs. (9), (13) and (14) are satisfied. 
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B. Time transformations and expansion 

The densities e± in Eq. (20) are given in terms of co- 
ordinate times t± . However, it is the time t measured by 
an observer on the shell which appears in the evolution 
equation (5). Hence we need to find the transformation 
from t± to t. Following Ref [9], we recall that metric 
continuity (1) implies 2 







\Hpp 


~dr~ 


s 


\ A 



1± 



A 2 - A (p 2 + H 2 p 2 ) 



,(25) 



in which H is again the Hubble "constant" , A = nep 2 — 1, 
and the expression within braces must be estimated on 
the two sides of the shell 3 . Now, the above two equations 
should be solved along with Eq. (5), which makes it clear 
why it is impossible to obtain general analytic solutions. 

An important result can be obtained by considering 
the time when the bubble is at rest, that is t± — to± and 
t = t , with p(tq) = and H = H(t ), namely 



dt± 
~dr~ 



±4 



Hi 



o± Po 



E,0 



*0± 



(26) 



From Eq. (16) we see that A = H 2 p 2 — 1, therefore 



dt± 
dT 



±1 



±1 



o± Pa 



±1 



V-Ao± ' 



(27) 



with the signs in the numerator simply reflecting the di- 
rections t± flow relative to r. It is clear that real solutions 
to Eq. (27) exist only if 



A ± < . 



(28) 



Remarkably, this is the same as stating that the en- 
ergy density inside the radius p = po must not gen- 
erate a black hole, as one can easily check by con- 
sidering the Schwarzschild radius r$ = 2 Gn M with 
M = (4 7r/3) £ /Oq. This is manifest when considering 
eo- inside the bubble, but it must also hold for the en- 
ergy density eo+ outside the bubble. For the outer region, 
this means the bubble must lie inside the Hubble radius, 
po < Hq 1 . Putting together these conditions tells us 
that the temporal coordinates are properly transformed 
only within a causal region of the spacetime. 

In order to study how the bubble grows after nucle- 
ation, we can expand t = t(r) for short times about to± 



2 There is a typo in Eq. (BIO) of Ref [9]: p 2 is missing in the last 
term in the square root. 

3 Note the sign ambiguity ± in front of the square root just reflects 
the double root of a second degree equation and is not associated 
with the interior or exterior regions. 



and To- Further, we want all times directed the same 
way, so we choose the + sign in the above expression and 
obtain, to linear order, 



t± — to± + 



TO 



(29) 



where t ± are integration constants. 

Unfortunately, a first order expansion is not sufficient 
to study the evolution of the bubble radius. Since po = 0, 
we need at least second order terms in r to get significant 
results, which makes all expressions very cumbersome. 
We shall therefore just consider a few specific cases, gen- 
eralizing the exact result (8) for a shell of constant sur- 
face energy in vacuum. For such cases, our perturbativc 
approach will yield exact conditions for the bubble's ex- 
istence, which we see as a clear advantage with respect to 
purely numerical solutions. Other advantages would be 
that having analytical expressions is a necessary ingre- 
dient for quantum mechanical (or semiclassical) studies 
of these systems. Moreover, adapting our procedure to 
all possible combinations of fluids in fl±, and for more 
general shell surface density, should be rather straight- 
forward. 



III. RADIATION BUBBLE INSIDE DUST 

The main idea in our approach stems from the observa- 
tion that the (three) fundamental (first order differential) 
equations (5) and (25) contain six functions of the proper 
time r: the shell radius p, its surface density a, the two 
times t± and the two Hubble functions H±. Once we 
choose the matter content inside fi± and on the shell, 
the Hubble functions and surface density are uniquely 
fixed and we are left with the three unknowns p and t± 
(and a set of constraints for the initial conditions). To 
determine these unknowns, we find it convenient to for- 
mally expand the shell radius p and Hubble functions H± 
for short (proper) time "after the nucleation of the bub- 
ble" (when = 0), and solve Eqs. (5) and (25) order by 
order. 

Since expressions rapidly become involved, and a gen- 
eral treatment for all combinations of matter content in 
0± and shell surface density would be hardly readable, 
we shall only consider the specific case of a radiation bub- 
ble (w- = 1/3) inside a region filled with dust (w + = 0). 
We further assume the shell's surface density 



ct(t) = cr > 



(30) 



is constant and positive. Since one would expect the 
shell's density decreases as the shell's surface grows, this 
assumption might appear rather strong. However, it is 
the simplest way to ensure the junction remains of the 
"black hole" type, and a more thorough discussion of 
this point can be found in Ref. [9]. In order to keep 
the presentation uncluttered, we also set k = 1 from 
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now on and regard all quantities as dimcnsionless (tan- 
tamount to assuming they are rescaled by suitable pow- 
ers of k = £p/Mp). This means that densities will be 
measured in Planck units, that is e = 1 corresponds to 



the Planck density e P = M P /i 



and 



1 to 



cr P = M P /£ 



Likewise, p = 1 is the Planck length 



l-p. We also express the shell surface density and radia- 



tion energy density as fractions of eQ USt > 0, 



rad _ dust 
e — e x i 



(TO 



, /.-dust 

- V ° 



y 



(31) 



with < x < 1 and y > 0. 

It is natural to choose To = 0, and then proceed to ana- 
lyze Eqs. (5) and (25) by formally expanding all relevant 
time-dependent functions in powers of r — To = r: 

Step 1) since p = [sec Eq. (9)], we can formally 
write the bubble radius as 



p(t)= Pq + ^p t 2 + O(t 3 ) 



(32) 



where po and po are parameters to be determined. In 
particular, from Eqs. (6) and (10), we obtain the (not 
yet final) expression 



Po 



3 cto 



\j (eo+ + eZ + 9 « erg/ 4 ) 2 - 4 e o- 

3 ( e dust)-l/2 y 



(33) 



Ax 



which only depends on eo± and ao- More precisely, ejj ust 
sets the overall scale of the shell radius and the fractions 
x and y defined in Eq. (31) the detailed form. 

We next obtain t± — t±(r) by solving Eq. (25). How- 
ever, for this purpose we need the Hubble parameters H± 
as functions of r, whereas they explicitly depend on t±: 

Step 2) we replace H in Eq. (25) with the formal 
expansion 



H = Ho+Hot + 0(t 2 ) , 



(34) 



where Ho and Ho are unknown constant quantities to be 
determined by consistency. By expanding the right hand 
side of Eq. (25) to first order in t and then integrating, 
we obtain t to second order in t, 



t± — to± + 



7=A 



+ 



o± 



PqHq± 
2 A, 



o± 



Po 



PoH 



0± 



! ,(35) 



o± 



where po must now be understood as the expression given 
in Eq. (33) and t ± are integration constants we can set 
to zero without loss of generality. In fact, let us assume 
we are at rest in an expanding (or contracting) universe, 
corresponding to the old exterior phase (with parameters 
e 0+ and H +), and measure a time t' + from the "Big 
Bang" (or beginning of collapse) of this exterior universe 
[a+(t' + = 0) = or a+(t' + = 0) > a(t' + ) for t' + > 0, 



respectively]. If, for instance, a new phase bubble arises 
at rest at the instant t' + = t' 0+1 we can define t + = t' + — 
t' 0+ , so that the bubble is created at t+ = 0, and also 
set t - = 0, since an "inner time" is meaningless before 
any "inner part" exists. For different pictures, similar 
arguments can likewise be formulated. 

Step 3) From Eq. (21) and (22), we choose an expand- 
ing radiation interior and contracting or expanding dust 
exterior, 



= 4 ad (M , 



a+ = 4 t ust (t+) 



(36) 



set a ± = 1 and express t± according to Eq. (35). In 
so doing, a± and da±/dt± become explicit functions of 
t containing p , Ho± and Ho±- For consistency with 
Eq. (34), we must therefore require 



H(r) = 



1 da 
a dt 



H(i 



(37) 



with H- = H- > and H + = H+ < for collapsing 
dust, or H + = H+ > for expanding dust. 

Step 4) To zero order in r, Eq. (37) gives rise to a 
first-order equation for Ho, 



Ho 



_L_ da 

a dt 



T = 



da 
d* 



Ho ■ 



T = 



The solutions are uniquely given by 



H 



n 



(38) 



(39) 



in which the f and + sign (respectively 4- and — sign) refer 
to expanding (contracting) solutions, i.e. solutions with 
increasing (decreasing) scale factor, as before 4 . Note 
this result also follows directly from the Friedmann equa- 
tion (16) for r = t = 0. 

To first order in t, one analogously obtains 



n — - 



ne 



Po £ o 



(40) 



with n = 2 for radiation and n = 3/2 for dust, and 
Po must again be understood as the expression given in 
Eq. (33). 

Step 5) Replace po from Eq. (33) and the chosen 
combination of Hubble parameters (39) inside Bq, which 
must then satisfy Eq. (13). This equation will only con- 
tain eo_ = e ad , eo+ = eQ USt and ao, so that it can be 
used to determine ao = cro( e o-, e o+)- ln particular, in- 
troducing the fractions in Eq. (31), we obtain 



Bo 



^dust 



bo{x,y) 



(41) 



4 Note that, for example, the Hubble parameter for the expanding 
interior phase will carry a second subscript sign and will then be 
denoted as H_ , where the subscript — indicates the inner region 
and the apex f stands for expansion. 
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from which it appears that the dust energy density just 
sets the overall scale like in Eq. (33). For any given values 
of e dust , the shell surface density is instead determined by 
the radiation energy density according to 



bo(x,y) = , 



(42) 



which, for the cases of interest, is a fourth-order algebraic 
equation for y. Analytic solutions can be found (in suit- 
able ranges of x), which we denote as y = y(x), so that 
the allowed surface densities are given by 



(To 



-dust 



y 



(43) 



Step 6) Replace the above surface density Co into the 
initial radius (33) and obtain its final expression, 



Po 



3 ( e duBt)-i/2y 



(44) 



which can then be used to determine the final forms of 
Hq± and the scale factors a± to first order in r. 

Step 7) One must now check that cto and po satisfy all 
of the initial constraints and lead to valid time transfor- 
mations (25), at least for some values of e ad and eQ USt . If 
not all of these conditions are met, one must conclude the 
corresponding physical system may not exist. Moreover, 
we note the condition (28) requires eo < ep in order to 



have a (semi)classical bubble with po 



> 



In the follow- 



ing Section, one should therefore consider only dust and 
radiation energy densities eo <C ep and look at the limit- 
ing case eo — ep as a glimpse into the quantum gravity 
regime. 

If a consistent solution for tro and po exists, one can 
proceed to determine higher orders terms (in r). How- 
ever, due to the increasing degree of complexity of the 
resulting expressions, we shall not go any further here. 
We instead present our findings for the two cases of in- 
terest separately. 



A. Collapsing dust 

This case is defined by choosing the scale factors 

a- = 4 ad (<_) , a+ = a dust (i+) , (45) 

and proceeding as described above. We can then prove 
that this case does not admit solutions, in general, by 
simply analyzing the constraint (42), 

b oc 16x 2 (s + Ax 3/2 ^ +4x 3/2 (4-9y 2 )(4-8x-9y 2 ) 
+3(4 + 9y 2 )(4-8a; + 9y 2 ) =0 , (46) 
which admits the four solutions 



2 / 2z 3 / 4 ±iV3 



(47) 



FIG. 1: Plot of cr /\/eo ust = V+-( x ) (magnified by a fac- 
tor of 10 for convenience, dashed line) and corresponding 
Po/(eo us T 1/2 (solid line) in the range (53). 




FIG. 2: Plot of to+ (solid line) and to- (dashed line) for 
y — (x) in the range (53). 



However, for x ^ 1, all the y±± are complex and a com- 
plex surface density is obviously unphysical. One is then 
apparently left with the only trivial case x = 1, corre- 
sponding to y = and 



- o = (£ d U st r l/ 2l 



'3 + 4a; 3 / 2 
3 + 4x 1 /2 



1 . 



(48) 



This case however does not satisfy all the required con- 
straints. For example, Eq. (14) for Co = yields 



^dust 



>e^ d = a;e: 



dust 



(49) 



which clearly contradicts x = 1. Correspondingly, the 
time transformations (25) are not well-defined, because 
to+ = i+(r = 0;x) is complex for < x < 1 and both 
to± diverge for x — > 1. 

The overall conclusion is thus that expanding radiation 
bubbles with a turning point of minimum radius cannot 
be matched with a collapsing dust exterior. 



B. Expanding dust 



This case is defined by choosing the scale factors 



4 ad (t_) 



_,dust 



(M 



(50) 
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FIG. 3: Plot of po for y = y+-(x) with eg ust = e P /10 (dashed 
line) and eo ust = £ p (solid line) in the range (53) . Only values 
above £p = 1 represent acceptable semiclassical radii. 




10 x 



The only solution which appears consistent is therefore 



„ /.-dust 



y+- 



x) 



2x3/4 + ^ ' 



(55) 



with x again in the range (53). This expression yields 
a vanishing surface density for the limiting values x — > 
1 and x — > x m ; n (see Fig. 1) and further satisfies the 
condition (14), 



■vl- 



V3(l-x) 

2x 3 / 4 + 



> . 



(56) 



The corresponding initial bubble radius is an increasing 
function of x (see Fig.l), 



Po = (C St ) 




dust \ -1/2 



< (C st ) 



(57) 



with p (z -> 1) = (e dust )~ 1/2 - Further, the products 



e d ast p 2 <l and e^p^l 



rad -2 



(58) 



FIG. 4: Plot of C /M$ ust for a > and y = y+-{x) in the 
range (53). 



The crucial task is again to solve the constraint in 
Eq. (42), namely 

b oc 16 x 2 (3 - 4 x 3/2 ) - 4 x 3/2 (4-9 y 2 )(4 - 8 x - 9 y 2 ) 



-3(4 + 9y 2 ) (4-8x + 9y 
admitting the four solutions 

y±± =±l\ (i-x) 



2x3/4 + ^/3 
2x3/4^^3 



which are real for 



2/3 



< x < 1 . 



(51) 



(52) 



(53) 



We discard the negative solutions y_± associated to neg- 
ative surface densities and just analyze the positive cases 
y+±. It is then easy to see that y++ leads to a surface 
density that diverges for x — > x m ; n , and is always too 
large to satisfy the condition (14), since 



V3(x- 1) 



4 y++ 2x 3 /4 -VI 



<0, 



(54) 



in the range (53). In the limit for x — >• 1, y ++ — >• 0, 
however, Eq. (14) is still violated in the strict sense and 
one can in fact show that io+ diverges. 



for x m j n < x < 1, as required by the condition (28). 
In fact the initial time derivatives to± are well defined 
in this range (see Fig. 2) and only diverge for x — > 1. 
Note the above initial radius can be larger than £p only 
if eQ USt < ep and for sufficiently large x, since po — > for 
x — > x m ; n (sec, for example, Fig. 3). 

Finally, let us check if one can use the process of bubble 
nucleation to describe a phase transition from dust to 
radiation for the matter inside the sphere of radius po> 
accompanied by the creation of a layer of non-vanishing 
surface density ctq. From Eqs. (55) and (57), one has 



C = M. 



dust 



M rad 



Ek <0 



(59) 



which means the dust energy inside the sphere of radius 
po at time of bubble formation, M dust = (4tt/3) p\ e dust , 
is not sufficient to produce the radiation energy M™ d = 
(4 7r/3) po eQ ad and surface energy Eq — 4 tt <tq p\. An ex- 
tra source is thus needed to provide the energy —Co > 0. 
The reverse process of collapsing radiation reaching a 
minimum size p — po and then turning into collapsing 
dust would instead be energetically favored, with the 
amount of energy —Co now being released. Of course, 
in order to support this kind of argument, the extra con- 
tribution should be a small perturbation on the given 
background, 



|Co 



< M dust 



(60) 



since it was not included in the dynamical equations. 
From Fig. (4), we expect this is indeed a very good ap- 
proximation since < —Co < 0.06M dust . 
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IV. CONCLUSIONS AND OUTLOOK 

In this paper, we have analyzed bubbles of radiation 
whose time- like surface starts to expand inside collapsing 
or expanding dust with vanishing initial rate, and with 
the further (simplifying) assumptions that the bubble's 
surface density is constant and positive, and both interior 
and exterior are spatially flat. These bubbles generalize 
the simplest self-gravitating case of a shell with constant 
surface density expanding in vacuum, for which the ex- 
act trajectories are known [1, 3]. These generalizations 
arc of potential interest both for the physics of the early 
universe and the description of astrophysical processes. 
However, although the general formalism was already de- 
veloped a long time ago [1], and the dynamics are ruled 
by apparently simple equations [9], finding explicit solu- 
tions is not straightforward. 

By developing an approach to obtain analytical ex- 
pressions for the evolution of the bubble radius in the 
shell's proper time, p = p(r) with p(r = 0) = 0, we 
determined the conditions which allow for the existence 
of such configurations. Although our approach is per- 
turbative (with an expansion for short times after nucle- 
ation), the conditions for the bubble's existence are exact, 
which is a clear advantage with respect to purely numer- 
ical solutions. We then found that expanding radiation 
bubbles of constant surface density may not be matched 
to a collapsing dust exterior. More precisely, we found 
that inside collapsing dust there may not be a bubble of 
radiation whose surface ever reaches vanishing speed of 
expansion at finite radius. Bubbles whose radius admits 
a turning point are instead allowed inside an expanding 
dust-dominated universe. They further can be used to 
model a phase transition from radiation to dust if an ex- 
ternal source provides part of the the energy required to 
build the shell, or the converse process with release of 
energy (albeit, of an amount small enough to leave the 
background configuration unaffected). 

Let us clarify this point about energy conservation. 
The fundamental Eqs. (5) and (25) are just a different 
form of the junction equations (1) and (2) which, in turn, 
follow from the Einstein equations. Conservation of the 
energy-momentum in a given system is therefore guaran- 
teed. However, when we use bubbles to model a phase 
transition, we are considering the possibility that a region 
of space filled with dust be replaced by radiation enclosed 
inside an expanding shell, or the reverse process. Techni- 
cally, we are therefore considering two different systems: 
one with dust and one with a bubble of radiation within 
a shell of positive surface density whose radius evolves 
along a trajectory with a turning point (zero speed at 
finite minimum radius) . The total energy in the two con- 
figurations differ by the amount Co defined in Eq. (59), 
and a (quantum) transition between them would there- 



fore violate energy conservation and be suppressed in the 
scmiclassical regime. By looking at Fig. 4, we however 
see that |Co| <C Mq us *. One may thus argue the unspec- 
ified matter contribution carrying the energy |Co| should 
be well approximated as a perturbation with respect to 
the dust and radiation, with its backreaction on the cho- 
sen configuration consistently negligible. If so, one can 
further speculate if the extra energy required to nucle- 
ate radiation could be provided by pressure in the initial 
cloud or by the decay of a region of false vacuum (with 
vacuum energy or cosmological constant A + ) to true vac- 
uum (with cosmological constant A_ < A + ), like in the 
seminal Refs. [4]. 

The fact that no consistent solution was found inside 
collapsing dust does not mean that expanding radiation 
bubbles may not be produced at all in this context, which 
includes, for example, the collapsing core of a supernova 
or other astrophysical processes leading to black hole for- 
mation. In fact, the situation might change if one, for in- 
stance (and more realistically), includes matter pressure 
or a radius-dependent surface density, a = o{p). This 
observation thus brings us to briefly comment on the pos- 
sible generalizations and extensions of the present work, 
which include the just mentioned non-constant cr, as well 
as different combinations of matter inside and outside 
the shell, and cosmological constant (s) A±. Moreover, 
one might like to consider the vacuum inside the shell 
and radiation outside (with or without A±) and apply 
the corresponding results to the thick shell model previ- 
ously studied in Refs. [12]. 

Finally, our analysis is entirely based on classical Gen- 
eral Relativity and no attempt was made to compute the 
quantum mechanical "tunneling" probability for radia- 
tion bubbles to come into existence (or convert to dust). 
Such an analysis requires the (effective Euclidean) action 
to be integrated along the (classically forbidden) trajec- 
tory for the bubble radius p to go from to po [5] , whose 
construction is clearly no easy task, given the classical 
trajectories are so difficult to determine. Nonetheless, 
another advantage of our approach is that it provides 
analytical (albeit perturbative) expressions, which is a 
property one needs for any quantum mechanical (or semi- 
classical) studies of these systems. Of course, energy den- 
sities above the Planck scale would not be meaningful in 
this context, since one then has no guarantee the dy- 
namical equations derived from General Relativity can 
be trusted at all. 
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